Electron density fluctuations accelerate the branching of streamer discharges in air 
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Branching is an essential element of streamer discharge dynamics but today it is understood only qualitatively. 
The variability and irregularity observed in branched streamer trees suggest that stochastic terms are relevant for 
the description of streamer branching. We here consider electron density fluctuations due to the discrete particle 
number as a source of stochasticity in positive streamers in air at standard temperature and pressure. We derive 
a quantitative estimate for the branching distance that agrees within a factor of 2 with experimental values. As 
branching without noise would occur later, if at all, we conclude that stochastic particle noise is relevant for 
streamer branching in air at atmospheric pressure. 



Streamers are filamentary electrical discharges that propa- 
gate through a non-conducting medium when it is suddenly 
exposed to a high electric field ||T|. As first stages of electri- 
cal breakdown |2|, they are found in nature preceding a light- 
ning stroke and as the building blocks of upper- atmospheric 
discharges above thunderstorms OJjS]. Due to their efficient 
production of chemical radicals f6l, streamers are used in in- 
dustry for gas cleaning and sterilization. 

Streamer discharges frequently form irregular trees with 
many branches both in laboratory, where branch length and 
branching angles in air have been measured recently Q (8], 
and in upper- atmospheric discharges |9|. However, our 
present understanding of streamer branching is only qualita- 
tive. Streamer trees in the laboratory are highly random and ir- 
regular but simulations have shown that streamers can branch 
even in fully deterministic density models lilOH13l . Model re- 
duction and analytical theory explained branching as a Lapla- 
cian instability that develops when the space charge layer 
around the streamer head is much thinner than the streamer 
radius il4j . Nevertheless, the analysis of a reduced moving- 
boundary streamer model suggests that streamer heads are lin- 
early stable ifTSl [T6l even if the stabilizing effects of electron 
diffusion and photo-ionization are neglected. Streamer fingers 
are therefore similar to laminar pipe flow: a finite perturbation 
is required to let streamers branch or to make the pipe flow 
turbulent. 

A small, but finite perturbation to trigger branching can be 
due to the initial condition. But it can also be created by 
fluctuations during the evolution. Such fluctuations are nat- 
urally included in Monte Carlo models |17-19| that follow 
the single electron motion; they model the stochastic distribu- 
tion of electron positions and energies. Those models are mi- 
croscopically very accurate at the cost of high computational 
demands. Even state-of-the art spatially hybrid codes 1 19 1 are 
still limited to quite short streamers. 

Here we introduce a new computational model and present 
its quantitative predictions for positive streamers in air at stan- 
dard temperature and pressure. The model is a spatially ex- 
tended stochastic model |20| that assumes an electron energy 
distribution determined by the local electric field but accounts 
for density fluctuations due to the discrete number of elec- 



trons. Since particles are not individually tracked, the required 
memory and computations are roughly independent of the 
number of active particles, in contrast to Monte Carlo meth- 
ods. We study the evolution towards branching with realistic 
density fluctuations in full three dimensions and, extrapolat- 
ing, obtain branching ratios that are consistent with the exper- 
iments |7, 8|. This suggests that streamer branching in air at 
atmospheric pressure is triggered by electron density fluctua- 
tions. 

Model.- The most relevant microscopic processes in 
streamers are two-body reactions between free electrons and 
neutral gas molecules. Therefore most quantities of a streamer 
discharge scale with the neutral gas density in a definite man- 
ner called Townsend scaling [5 J. Following 121] we define a 
typical length for streamers in air Iq a; 2.3 yum ■ (Nq/N), a typ- 
ical electric field 2 ■ 10^ V/cm ■ (Nq/N) a typical time 
fo = 3 ■ 10"'^ s ■ (No/N) and a typical density of charge car- 
riers no = £{)E()/el{) x 4.7 ■ lO^'^cm"-' ■ (No/N)^, where e is 
the elementary charge, is the molecule number density of 
air and A^o is the number density at standard pressure and tem- 
perature, used as an arbitrary reference. We refer to Q for a 
physical interpretation of these quantities and scaling laws. 

Using these typical magnitudes one can build a dimension- 
less streamer model. We consider here a minimal model for 
air 1211 l22j that also includes photo-ionization |231. We are 
interested on the dynamics of the streamer head, where im- 
pact ionization strongly dominates over electron removal by 
dissociative attachment and the latter can be safely neglected. 
The governing equations are therefore 

d,o- = DVV + V ■ (crE) + S^"''> +s'^""P""\ (1) 
d,p ^ S'-"''^ + S^'"""'"\ (2) 
V^(f)^cr-p, E = -V0. (3) 

Here cr and p are the (dimensionless) electron and ion densi- 
ties, S "^''''^ and S (""'""^'') are, respectively, the photo-ionization 
and impact ionization sources of electron-ion pairs, (p is the 
electrostatic potential E the electric field, and D is a diffusion 
coefficient, taken as D = 0.1. 

There are several corrections to Townsend scaling fS) . One 
arises from collisional quenching of photo-ionization |12|, 
here included in 5*'^''', which contains an explicit dependence 
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FIG. 1. Scheme of an elementary time-step in the lattice model. 
Each electron is represented here by a dot; note however that we 
do not keep track of every individual particle: only the number of 
particles in each cell. During the time-step from t to t + At we move 
Ujj electrons from cell C, to Cj. Meanwhile, r, electrons arise in Cj 
from impact ionization and photo-ionization. 



on the air density. A second correction, on which we fo- 
cus here, arises from the finite number of particles and the 
stochastic nature of microscopic processes; this cannot be 
expressed in equations ([T])-(|3]l, because they only contain 
macroscopic quantities. Rather we will take a discretization 
of ([T])-([3]l and convert it into a spatially extended stochastic 
model|20|. 

Method.- Since the importance of stochastic noise depends 
on the particle number density our first step is to derive, from 
the magnitudes described above, a dimensionless parameter 
for the typical number of charge carriers contained in a typical 
volume. We define g = «o/q = EaSol^/e. This is the number of 
elementary charges that has to sit in each area of an infinite 
charged plane to create a jump Eq in the electric field. In air 
g ~ 5700 ■ {N[)/Ny, at an altitude of 70 km, typical for sprites, 
g.lO^ 

The relative amplitude of the statistical fluctuations of a 
number g of particles is g^'^^. The limit of negligible fluctua- 



tions is <K 1. In air at atmospheric pressure, g 



-1/2 



0.01 



and stochastic noise is a relatively small correction on the fluid 
description. However, as we will see, due to the strongly non- 
linear nature of streamer discharges, such small fluctuations 
can be amplified by strong electric fields and alter signifi- 
cantly the propagation of a streamer For sprites, g"''^ is much 
smaller, about 10 

Now let us take a spatial discretization of ([T]i-([3]i; the sim- 
ulation domain (see Fig[T]l is divided into cells C,, each with 
a dimensionless volume u,- = Vj/Iq (V, is the dimensional vol- 
ume). If we are given the dimensionless densities in each cell, 
cr, and Pi, we know, from our discretization, how to calculate 
the left hand sides of ([T|-([3]l. In particular, we can calculate 
the source term S',- = Sf^^^ _^ ^(mpaa) ^.j^^ from cell C,- 
to each neighboring cell Cj, which we denote Fij. 

But instead of densities one has a discrete number of elec- 
trons and ions in each cell, Nf* and A^*'^; the dimension- 
less densities are therefore cr,- = Nf^ jnoVi - Nf^ Igvj, p, - 
N^'^ /gvi. We may now use these to calculate the source terms 
and the fluxes. The problem is now how to update the number 
of particles as time evolves. To simplify the description, let us 



first describe a time-stepping that will tend to a forward Euler 
discretization with time step Af. 

We look first at the source terms. Let r, be the number of 
electron-ions pairs created in C, during a time At. This quan- 
tity follows a Poisson distribution with average A, - gVjS ,Af . 
So in the simulation we draw for each / a sample from the 
Poisson probability distribution pir,) = /l^'e"*/r,!. 

The transport terms are slightly more complicated because 
the number of electrons must be conserved. We interpret the 
gviFjjAt as the average number of electrons that flows from 
C, to Cj during At (note that Fij + F ji). Thus the prob- 
ability for an electron in C, at t to end in Cj at f H- Af is 
Pij = NfFijAt/gvi for / ^ j. But since the electron must 
end somewhere, pa - I - Yjj^i Pij P24l. We can use these ptj 
to obtain the number of electrons moving from C, to Cj, that 
we denote u/j. The probability distribution for the electrons 
exiting C,- is the multinomial distribution with Yjj Uij - Nf^. 

Now we have all the ingredients to update the particle num- 
bers as 



Nfit + At) = Nf(t) + n + Y, - «o) 



(4) 



As mentioned above, this scheme tends to an explicit Euler 
time discretization of the fluid equations as g^^^^ 0. But 
we can also design a two-step time updating that tends to sec- 
ond order Runge-Kutta ll25l if we (a) use the particle numbers 
at t to calculate 5,(f) and Fij(t), (b) perform a half-step using 
Af/2 to update the particle numbers, (c) use the new parti- 
cle numbers to obtain S i{t + At/2) and S i(t + At/2) (c) define 
S'.{t) = (Si{t)+Si(t+At))/2, F'.jit) = {F^{t)+F^(t+At))/2 and 
use them to perform a step Af. Note that in principle step (b) 
could also implement a standard continuum step, since noth- 
ing forces us to preserve an integer number of particles in that 
intermediate step; we opted however to use the same stochas- 
tic step in both stages of the algorithm. 

We have not yet mentioned the spatial discretization to cal- 
culate Fij. The reason is that the scheme is flexible on that. 
We used here the scheme described in 1251 . This is a flux- 
limited, nonlinear discretization schemes and it poses an addi- 
tional difficulty: it sometimes leads to negative Fij which can- 
not be interpreted as a probability. In that case, we rearrange 
the fluxes by letting Fij — > F'^j, with F^ - 0, F'j. = Fji - Fij 
until no negative fluxes remain. 

Results.- We now use the adaptive grid refinement and the 
fluxes and reaction terms from |25 l and the three-dimensional 
cylindrical mesh of Ref. [26J. As a first application, let us 
analyze the initiation of breakdown in a small plane-to-plane 
geometry with a potential difference of 16 kV between two 
electrodes separated by 2 mm of air. The simulated volume 
is discretized into cells Ar = Az = 8 pm, AO - 2n/6A. As 
initial condition, we set a neutral hemispherical gaussian seed 
at z = r = (positive electrode) containing ~ 6 ■ 10^ electrons. 

Figure |2] shows the evolution of a cross-section of the elec- 
tron densities up to 1.35 ns. In that short time-span, a multi- 
tude of avalanches seeded by photo-ionization has developed. 
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FIG. 3. Highest deviation of tlie absolute value of the electric field 
from its azimuthal average at three snapshots. The volume labelled 
A corresponds to the evolution of the streamer up to ? = 13.5 ns, with 
noise switched off. When we switch on noise, the streamer body 
(B) develops low-amplitude, small-scale fluctuations. The streamer 
surface close to the tip (C) has fluctuations with a higher amplitude 
and an auto-correlation length of about one tenth of a millimeter. 
The amplitude grows but in this simulation it remains small small 
compared with values of some hundreds of kV/cm in that area. 



FIG. 2. Snapshots of the cross section of the electron density in the 
y = plane in a plane-plane, 2 mm gap with a potential difference 
of 16 kV. Multiple avalanches start out of seeds produced by photo- 
ionization; the electric breakdown extends to the complete volume 
and no actual streamer is initiated. A movie of this simulation is 
available in the auxiliary material. 

A very similar evolution is observed in the Monte Carlo sim- 
ulations of I.19J; both results show that in strong electric fields 
noise weakens or even prevents the formation of streamers. 

The situation changes however for a streamer in a point- 
plane geometry. Then only in a small volume is the electric 
field high enough to produce significant ionization. 

We run a simulation with point-plane electrodes imple- 
mented as in [27 1 . The needle is 2 mm long and its tip is sep- 
arated from the plate by 7.2 mm. As initial condition we set a 
semi-spherical neutral gaussian ionization seed at the needle 
tip with a radius of 73.6 jim and a peak ionization density of 
4.7 ■ 10'^ cm"^. The needle has a positive potential of 10.5 kV 
relative to the lower electrode. We used an adptive refinement 
strategy [25] with coarsest grid Az = Ar = 40 jum and finest 
grid Az = Ar - 2.5 jum |[28 l. Since full 3d simulations are 
too demanding, we chose to run the simulation with cylin- 
drical symmetry and g^^^^ = up to f = 13.5 ns. Then we 
remove the constraint of cylindrical symmetry and introduce 
stochastic noise at the level expected at atmospheric pressure 
ig = 5700). The continuous density at each cell is then in- 
terpreted as an average and discrete numbers of particles are 
obtained by drawing random samples from a Poisson distribu- 
tion with this average. 

To represent the evolution of this 3d simulation that de- 
viates only slightly from perfect cylindical symmetry, let us 
consider the average of some quantity around the azimuthal 



angle, {u)g - J^'^ uij, z, 6)d6. The deviation from symme- 
try can then be defined as 6u - max0(u-{u)g). Figure|3]shows 
<5|£'| at three instants of time after noise is introduced into the 
simulation. 

This simulation is highly demanding: the short streamer 
evolution represented in Fig.|3]took about 5 weeks using two 
dual-core 3 GHz AMD Opteron processors. We could not run 
the simulation long enough to observe actual branching; nev- 
ertheless, we can use the present result for a first quantitative 
estimation of the time needed to branch. 

Let us look at the Fourier transform of the electron 
density cr along the azimuthal coordinate 6, &{r,z,k) - 

Jq d9ne(r,z,0)e^'''^. We can define the "total spec- 
tral content" of mode k of the electron density as Wk = 

X-j-CO nCO 
^ dz Jq drr&{r,z,k). If the streamer is close to cylin- 
drical symmetry (i.e. it is far from a branching state), |Wo| » 
\Wk\ for all k > Q. Hence if we define A = 2i>o IW^lVl^ol^ 
we can postulate that the condition for branching is A ~ 1 . 

In Figure |4] the value of A during the simulated time frame 
is plotted as a black line. During the full evolution we find 
A <s 1, i.e. the streamer was at all times far from branching. 
However we can use A(f) to obtain a first estimation of the 
branching time. After a transient, A growths exponentially; 
extrapolating this growth to A ~ 1 we can roughly estimate 
the branching time as t branch ~ 21 ns. The streamer velocity 
is approximately u = 0.32mm/ns and hence after the intro- 
duction of noise the streamer would run for about 2.4 mm be- 
fore branching. The measurements in QE) give a ratio of the 
branching distance to the streamer diameter of about 12-15. 
This is relative to the radiative diameter, estimated to be about 
half of the electrodynamic diameter |29 1; with that estimation, 
our value is about 8. One must also take into account that (a) 
short branching distances are harder to measure and therefore 
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FIG. 4. Evolution of A as defined in the text. We show the evo- 
lution of three runs: in all runs the system was deterministic and 
cylindrically symmetrical up to time 13.5 ns. After this point we 
add stochasticity and allow deviations from cylindrical symmetry. In 
the curves marked as deterministic, we remove noise again at the 
time marked with the dots while keeping the unsymmetric perturba- 
tions that have evolved up to that time. Within this limited time- 
frame, A seems go through a transient, fast growth phase until it 
settles into an exponential growth. The dashed line fits this second 
phase. For the completely stochastic run the best fit is A = ae'''" with 
a = 5 • 10"^", T = 0.48 ns. Extrapolating, this predicts branching at 
/ rloga X 21 ns. 

the average in Q |8] may be slightly overestimated and (b) 
both in our model and in observations branching distance is 
random: we are comparing the result of only one simulation 
with an average over many measurements so a certain discrep- 
ancy seems natural. 

However, our estimation of a branching distance is close to 
the measured value, suggesting that noise resulting from the 
finite number of particles is relevant in streamer discharges. 
To measure the relevance of a persisting noise compared with 
the inherent growth of cylindrically perturbed modes we per- 
formed two simulations in which noise is removed after some 
time. These are shown in the two curves of Figure]?] labelled 
as deterministic (a) and (b). We see that noise always in- 
creases the growth rate of the deviations but also that even 
a relatively small amount of noise during a short time (a) is 
enough to trigger an instability that would eventually lead to 
streamer branching. Longer simulations must be performed 
to check that the evolution suggested in Fig.]4]continues until 
the streamer branches. 

Note that at lower pressure, such as those in the high layers 
of the atmosphere where sprite streamers are commonly ob- 
served, stochastic noise is weaker {g-^>^ * 10"'*). Although 
the ratio between diameter and branching distance has not 
yet been measured in sprites, they also branch frequently 19). 
However, in the upper atmosphere, other sources of stochas- 
ticity may be relevant, such as cosmic rays and atmospheric 
inhomogeneities. 
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